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Abstract : 


Extra terrestrial materials, certain materials which 
have their origin at the earth’s surface due to meteoritic 
impact, or under highly reducing conditions, such as in the 
case of basaltic flows in contact with coal beds or ser- 
pentenites, all contain Fe and FeNi phases with high mag- 
netization values and spherical shape. Normally, the demag- 
netizing field (Hp-NIg :where N is the demagnetizing factor 
and Is the saturation magnetization) is corrected for. 
However, in disperse systems, such as most natural materials, 
we do not make this correction, and instead the particle 
shape effects are analyzed in terms of the saturation fields, 
Hs=H- D =NIg, and the magnetization differences (Alg). Dis- 
crete size modes of superparamagnetic (SP), multidomain 
(MD), and single domain (SD) particles result in reduced 
coercive force (He), increase in the value Rr (=ratio of 
remanent coercive force, Hr, to H c ), and decrease in the 
value Rj (ratio of remanent magnetization, I R , to satura- 
tion magnetization, Ig). The main distinctions between 
the various natural materials can be made by this approach. 
Hysteresis loops for terrestrial basalts, Fe and Ni rods 
and spheres, chondrite meteorites, lunar samples, impac- 
tites, and chondritic fusion crust are presented. 
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Introduction : 

Natural materials, with the exception of iron 
meteorites and terrestrial ores are dispersed ferro- 
magnetics, i.e. grains of ferromagnetic material dis- 
persed in a nonmagnetic silicate matrix or other non- 
magnetic matrix. In extra terrestrial materials, im- 
pactites and certain terrestrial rocks, such as basaltic’ 
flows in contact with coal beds or trees and serpentenites, 
Pe and Pe alloy phases are found which have high magne- 
tization values and add a significant magnetic grain 
shape aspect to the sample. 

Collisions in outer space, impact at planetary sur- 
faces, atmospheric heating, and in situ oxidation of land- 
ed meteorites will affect the ferromagnetic assem- 
blage in these natural materials. Almost nothing is 
known about the magnetic property changes associated 
with such processes. The utilization of magnetic 
hysteresis analysis has great promise (a) in evaluating 
different meteorite classes, (b) in evaluating the black 
fusion crust on chondrite meteorites, (c) in evaluating 
the grain size and shape distribution in natural mater- 
ials, (d) in evaluating the magnetic coercivity distri- 
bution, and (e) in evaluating the thermal, oxidative, and 
shock alterations in natural materials. This paper will 
demonstrate the effective utilization of magnetization 
measurements in evaluating magnetic property changes 


- 2 - 



\ 

/ 

associated with thermo chemical and thermomechanical altera- 
tions of the magnetic phases in natural ferromagnetic 
materials. 

Measurements and Mineral Phases ; 

The principle reasons for making magnetic measure- 
ments are : (a) evaluation of type, shape, and size of 

ferromagnetic phases, (h) evaluation of chemical and 
mechanical alterations of the mineral phases, (c) estab- 
lishment of a thermomechanical classification of chon- 
drites, lunar rocks, and impactites , based on magnetic 
properties, and (d) evaluation of remanent magnetism and 
the possible information embodied in any component. The 
initial step in such a program is the characterization 
of the sample. Since FeNi particles are present in 
extra terrestrial samples the role of these phases 
must be evaluated. 

Magnetic hysteresis loops and magnetization curves 
were measured on a vibrating sample magnetometer located 
at the U.S. Steel Research Laboratory. Fields up to 16 
K.Oe. were used. Samples of various natural and man made 
materials were studied: polycrystalline steel spheres 

(SAE 52100), nickel wire (99.99% purity), shocked FeNi 
alloys (8% Ni - 164 kb shock), crossed nickel wires, two 
sphere configurations, chondrite fusion crust, carbonace- 
ous chondrites, chondrules from the Allende and Chainpur 
chondrites, 30 ordinary chondrites with varying degrees 
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of atmospheric degradation, shock, and shock heating. 
Published and unpublished data on the lunar rocks and 
fines are also considered. 

In the chondrite meteorites, ferromagnetic mineral 
phases are present which cover a Curie point range of 
/V500° C. The significant mineral phases influencing 
magnetization curves are FeNi. Various natural proces- 
ses, such as cooling, oxidation, high temperature altera- 
tions and shock induced deformation will alter the phase 
composition and structures and the effective particle 
size and shape. Previous work on chondritic meteorites 
by Stacey £t al. (1961) has shown that "Fe-jO^ 1 ' accounts 
for 100% of the moment in the Mokoia carbonaceous chon- 
drite, and kamacite with 5-6% nickel accounts for 80-90% 
of the moment in the chondrites Bjurbole, Mount Brown, 
Homestead, Barrata, and Cumberland Falls. The remaining 
20-10% of the moment is not accounted for. 

Impactites contain a wide range of ferromagnetic 
particle assemblages. Reference to Chao (1968) and El 
Goresy (1968) and numerous other papers in the volume 
Shock Metamorphi sm of Natural Materials (1968) (ed. French 
and Short) will explain the nature of the opaque mineral 
fraction in impactites. Lunar samples contain iron parti- 
cles ranging from 100 2. to millimeter sizes. Iron is 
present as spherules and other shapes. There are other 
phases which may be ferromagnetic but which have not 
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been studied in any detail . Meteorites, lunar samples, 
and impactites differ from normal terrestrial samples in 
the predominance of shape anisotropy due to iron and iron- 
nickel phases. Discrete mixtures of particles of varying 
coercivity and/or shape anisotropy cause reduced coercive 
force and large ratios (R h ) of remanent coercive force 
(Kneller and Luborsky, 1964; Meikeljohn, 1953) (H r ) to 
coercive force (H c ) in meteorites, lunar samples, and 
impactites. 

Particle Shape and Saturation Magnetization : 

A chondritic meteorite find is usually oxidized to 

a degree depending on local conditions and time spent on 

earth. If an Pe sphere is converted to FejC^ the shape 

field (Hs= 4 W Is) changes from Hs r* 7000 Oe to Hs* 

3 

a# 2000 Oe. This situation is illustrated in Figure 1 
where hypothetical magnetization curves for Fe and Fe 30 ^ 
are shown. The Hs fields, the slopes of the curves, the 
susceptibility, and the saturation magnetization values 
(Is) are indicated in the figure. Particle size and 
shape distributions will introduce curvature in the I vs. 
H curves for H < Hs. It is, however, quite clear that 
magnetic fields at least 7000 Oersted or greater are 
needed in studying the extraterrestrial samples. To com- 
pare iron bearing samples with terrestrial samples, the 
effect of the saturation magnetization differences (A Is) 
must be accounted for. 
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Experimental Results : 

Magnetic hysteresis loops and magnetization curves 
for artificial and natural materials are presented with 
examples chosen specifically to illustrate the effect 
of particle shape and composition of the ferromagnetic 
phases. 

Basalt - Sample CH21-001 is a fresh basalt dredged from 
the Mid Atlantic Ridge. The sample is typical of pillow 
basalts where a glass crust grades into an aphanitic to 
crystalline interior. Curve A ( Figure 2) corresponds to 
the glassy crust and curve B to the aphanitic interior. 
The slopes of the lines AA* and BB* are characteristic 
of the paramagnetic component in the samples whose slope 
can be evaluated by considering: 

= A X 10" 6 • BX + C(Y +Z) (1) 

where A, B, and C are empirical constants and X, Y, and Z 
are respectively the amounts of Fe^*, Fe^ + , and Mn^*. 

The amount of paramagnetic phase can then be evaluated by 

C = a + b • 10^ p£m 

where C is concentration, a and b are constants and 
is evaluated as in (l). The same argument holds for all 
materials. As can be seen in the figure the glassy crust 
has a significantly larger He value than the aphanitic 
interior which is to be expected since the magnetite 
grain size distribution varies according to the cooling 
rate, and the smaller mode is skewed to the glassy crust. 
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The total. amount of ferromagnetic in the aphanitic inter- 
ior is of course much greater than in the glassy crust, 
which is reflected in the amount of paramagnetic in the 
glassy crust. The overall chemical composition of the 
glassy crust and aphanitic interiors are identical. 

Nickel and Iron - Magnetization curves for small rods 
and spheres illustrate the shape fields (Hg), the influ- 
ence of the demagnetisation factor (N), and the influ- 
ence of saturation magnetization ( 13 )* Magnetization 
curves for a 5mm long Ni wire measured parallel (a) and 
perpendicular (b) to the rod axis are shown in Figure 3A . 
The slope of the magnetization curve is l/N and is the 
susceptibility in terms of the applied field and is con- 
stant until at H^=NIg saturation is achieved. When wall 
coercivity is considered (H w ) then the saturation field 
is Hg=NIg+H w . 

When two wires (each 5.5 mm long) are crossed the 
results are quite different as shown in Figure 3B for 
measurements along the direction (c) parallel to one 
axis and perpendicular to another, and direction (d) at 
45° to the perpendicular wires. Both measurements were 
made with the applied field in the plane containing the 
wire axis. Both demagnetizing fields are seen by the 
magnetometer, e.g., the results shown in Figure 5A . 

This simulation is important since it has been shown to 
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be similar to the fusion crust in certain chondrite 
meteorites (Wasilewski, 1972) where sulfide veins criss 
cross and/or Fe^O^ coexists with metal. 

The results for the iron sphere and needle are 
shown in Figure 4A . The influence of the shape field can 
be seen by comparing Figure 4A and Figure 3A . If Ni and 
Fe spheres were mixed together, a resultant curve simi- 
lar to the crossed Ni wire case (Figure 3B) would be 
measured since the shape fields for particles with con- 
stant N depend on Ig. 

Impactite - In the Monturaqui impactite which contains 
metal grains, as well as magnetic sulfides and oxides, 
the resultant hysteresis loop is much more complex than 
the cases in Figure 3 and 4A. In Figure 4B the hystere- 
sis loop for the impactite is shown. Because of the 
wide variation in Ig and N for the different metal, sul- 
fide, and oxide phases, continuous curvature in the loop 
is observed for all up to about 6000 Oe. where satura- 
tion is achieved. 

Particle Agglomeration - The value of N for a sphere is 
4W*/3 and for a needle of infinite length N is 2"W . If 
the spheres are noninteracting then N=47T/3 and the shape 
field is 4 W Ig/3. However if the spheres interact 
then the field Hg is AW Ig(l-f) f=packing fraction . 

If there is any alignment, then the direction of the ap- 
plied field is important. In Figure 5 are hysteresis 
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loops for a two sphere configuration where is applied 
parallel (A) and perpendicular (B) to the axis of the 
two sphere configuration. The resultant Hg fields are 

h s (a)<h s (s)<h s (b) 

where H S (S) is the shape field 4T/3 Ig of a single sphere. 
Chondrites — Two L6 chondrites (Wasilewski, 1972) Modoc 
and Holbrook were studied. Spherical particles and 
particles of irregular shape are present in both samples. 
Hysteresis loops for both 16 chondrites are shown in 
Figure 6Ao The Rjj value (ratio of H^ to H c ) is 14.6 
for Modoc based on Hc=75 Oe. and Hr= 1100 Oe. with Rj 
(ratio of remanent magnetization to saturation magneti- 
zation) =0.03. The high R H value is due to the presence 
of spherical shaped particles and in general to particles 
of high Nig values. A steel sphere (52100) has an R H 
value of 16. 

Taylor and Heymann (1969) have defined shock classes 
for chondrite meteorites based on metallographic studies. 
Hysteresis loops for chondrites which fit three of their 
shock classes are illustrated in Figure 6B Richardton 
(Class A), Holbrook (Class B), Calliham (Class D) . 

These studies are preliminary and though indications are 
that the shape fields, H s , change with increasing shock 
alteration more experimentation is required. 

Fusion Crust and Oxidation - Previous research into the 
magnetic properties of the black fusion crust on chondrite 
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meteorites (Wasilewski, 1972) indicates complete 
mineralogical changes and associated thermo chemical as 
well as thermoremanent magnetizations. The example 
shown in Figure 1, where Fe is oxidized to FejO^, 
applies. The R H and Rj ratios and the field Hg all 
vary in predictable ways. Any constriction in the chon- 
drite hysteresis loop is also eliminated in the fusion 
crust. 

Shown in Figure 7 are hysteresis loops for the Bowes- 
mont chondrite and the weathered Bowesmont fusion crust. 
Percentage wise there is a stronger paramagnetic compo- 
nent in the meteorite than in the fusion crust which is 
to be expected since Fe in the silicate phase is oxidized 
in the crust, rendering the silicate mass devoid of para- 
magnetic iron. 

Lunar Samples - Nearly all the lunar samples contain sup- 
erparamagnetic metal. This results in reduced H c , increase 
in the R^ value, and constriction in the hysteresis loops. 
The presence of large multidomain components also reduces 
H Q and results in an increase in the R^ value. Paramagne- 
tism due to iron in glass and silicate crystals is quite 
significant in the lunar samples. Shown in Figure 8 are 
hysteresis loops for lunar samples 10021, 10085, and 
10084 (T. Nagata and F. Schwerer private communication). 

The saturation field (Hg) is ^ 7000 Oe. for the fines 
(10084) , *** 4000 Oe. for the breccia (10085), and^ 2000 Oe. 
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for the breccia powder (10021). 

Remanence Ratio and Coercivity - The values of H c , H R 
and the ratios R^, Rj for some of the samples shown in 
Figures 2 to 8 are listed in Table I. It has been shown 
previously that a classification of natural ferromagnetic 
materials can be made by utilizing the magnetic hysteresis 
data (Wasilewski , 1972). 

Titanomagnetite (xFe 9 «TiO ^ -l-xFe ^O/| ) - For a given particle 
shape, N is the demagnetization factor, and Hg (=Nl0) de- 
pends on Ig„ For the titanomagnetites the Ig value ranges 
from Ig/* 0 for x =0.8 to Ig^480 emu/cc for x=0, at room 
temperature. For spheres, needles, plates, and ellip- 
soids the demagnetization factors are well established, 
but for natural materials the grain shapes are more severe. 
In particular, in ocean and continental basalts, meteori- 
tic fusion crust, and spherules of cosmic origin dendri- 
tic grain shapes are common. Rhodes et al. (1962) have 
shown that the H c of dendritic crystals will be smaller 
than the H c estimated from dimensions of an equivalent 
ellipsoid. It has been established that the Rj values 
increase and R^ values decrease as the grain size decreases 
for magnetic class I basalts. The influence of composi- 
tion and size of titanomagnetites on the hysteresis 
loops is demonstrated in the Figures 9A and 9B . Hystere- 
sis loop squareness depends both on grain size and com- 
position. 
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Shape Dispersions ; 

A bimodal mixture of spheres and needles can he 
distributed in an inert matrix in four ways ( Figure 10A ) . 
To a first approximation mixture (C) represents a chond- 
rite meteorite. Neglecting the size distribution, all 
chondrite meteorites consist of a heterogeneous mixture 
of spherical shaped grains and irregular shaped grains 
of metal. (In the carbonaceous chondrites this argu- 
ment does not hold because of the significance of the 
oxide phases..) Each of the cases (a) to (d) give mag- 
netization curves distinctly different from each other 
(Figure 10B) . 

Consider case (a) and (b) which are identical ex- 
cept for a 90° rotation with respect to the measur- 
ing field. The effect due to the spheres remains the 
same but in case (a) the field sees the hard rod direc- 
tion and in case (b) the easy direction. The measured 
curve is a composition curve for a simple two component 
assemblage ( Figure 11 ) . 

From a composition curve of this type, assuming a 
simple two component assemblage, a quantitative assess- 
ment of the distribution of components, component aniso- 
tropy, and the shape fields can be made. B and B* cor- 
respond to the field Hg (=NIg) associated with the 
needle and A and A' correspond to that for a sphere. 

The sharpness of the discontinuities at B and A depend 
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on the behavior of the two components in the region 
of approach to saturation and the specific shape dis- 
tributions. Curvature in the regions 0 to B and B to 
A correspond to the angular dependence of the needle 
component with respect to the field and to deviations 
from sphericity in the case of the spherical component. 

The slope beyond A is due to two causes: (a) apparent 

paramagnetism due to finite measuring fields, and (b) 
a true paramagnetic matrix component. In decomposing an 
unknown curve of the type shoen in Figure 11 the para- 
magnetism is evaluated (= Xh) as the dashed line, OC, 
which has the same slope as A !, AA ,,, . The linear extra- 
polation to the line H=0 at A" marks the saturation mag- 
netization of the system as Msp+Mn, where subscripts sp 
and n correspond to sphere and needle respectively. The 
extrapolation ABB" allows separation of the total moment 
into contributions from (a) spherical particles and (b) 
needles. If the grain shapes are randomly dispersed, then 
some distribution function must be introduce^ and the dis- 
cretness of the Hg field disappears. 

Discussion and Conclusions : 

The magnetic hysteresis loops for lunar fines and 
lunar breccia powder ( Figure 8) are quite different from 
each other. The Hg field for the fines is about 7.5 K.Oe. 
and about 2 K.Oe. for the breccia powder. This type of 
analysis does not enable us to evaluate the superparamag- 
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netism 


. The following can be said about the samples: 

A. Paramagnetism is greater in the breccia sample. 

B. Spherical particle shape is a much more signi- 
ficant ferromagnetic component in the fines. 

C. Constriction, i.e. reduced coercivity, is found 
in both samples and is due to mixed components. 

A more detailed analysis should provide much more useful 
information. 

Mixtures of discrete size modes of super paramagnetic 
single domain or multidomain particles, or the presence 
of the same size modes but different magnetization values 
will produce significant alterations to the hysteresis 
loop. The presence of superparamagnetism and large multi 
domain grains will reduce H and result in loop constric- 

V 

tion. The shape field or saturation field is defined as 
Hs=Nl5 and the presence of discrete shape modes can be 
easily detected. Likewise for the same shape but dif- 
ferent magnetization values, the difference Ig) in 
magnetization values can be detected (N<&Is)» Any thermo 
mechanical alteration which changes the magnetization 
values, the shape distribution, the size distribution, or 
induces structural transformations can be effectively 
studied by evaluating magnetization curves and hysteresis 
loops. 
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TABLE I 


Sample 

H c 

Sr 

— H 

Hi 

Note 

SAE 52100 

50 

800 

16 


steel sphere 

Monturaqui 

90 

225 

2.5 

.07 

impactite 

Modoc 

75 

1100 

14.6 

.03 

L6 chondrite 

Calliham 

165 

270 

1.6 

c— 

CM 

• 

metal shock 
class D 

10021 

24 

410 

17.1 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


- Simplified magnetization curve for an iron 
sphere and magnetite sphere. The shape 
field (H s =4/3 7T*Iq) and magnetization 
values are indicated. 

- Hysteresis loops for specimens from a sam- 
ple of quenched ocean basalt. As one proceeds 
inward from the glassy crust to the aphanitic 
and crystalline interior the paramagnetism, 
coercivity, and hysteresis change in systemmatic 
fashion. 

- Magnetization curves for nickel needles and 
crossed nickel needles. 

A. Magnetization curves parallel and per- 
pendicular to the axis of a 5 mm needle. 

B. Magnetization curves (a) along the axis of 
one needle and perpendicular to the other 
for the crossed needle set and (b) at an 
angle 45° to mutually perpendicular needles. 

- Hysteresis loops for iron needle, steel sphere, 
and Monturaqui impactite. 

A. Hysteresis loops for steel sphere and for 
the case parallel to the needle axis. 

B. Hysteresis loop for Monturaqui impactite. 

- Hysteresis loops for a two sphere configuration 
measured (a) parallel and (b) perpendicular to 
elongation. This demonstrates the effect of 
interaction on the shape field. Similar inter- 
acting configurations are found in lunar samples 
and chondrite meteorites. 

- Hysteresis loops for chondrite meteorites. 

A. Hysteresis loops for L6 chondrites , 

Modoc and Holbrook. 

B. Hysteresis loops for meteorites with pro- 
gressive shock metamorphism in the metal 
phases. A,B,D on curves refer to Taylor 
and Heymann (1969) metal shock classes. 

- Hysteresis loops for the fusion crust and 
meteorite body - Bowesmont chondrite. The 
ferromagnetic component for the crust is 
(FeNi )^04 and for the meteorite, FeNi alloy. 

This analysis reveals that oxidation of FeNi 
to (FeNi^O^ shifts the shape field Hs, in- 
creases Rj, and decreases Rpj. 
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Figure 8 - Magnetic hysteresis loops for lunar samples 
' (T. Nagata and F. Schwerer), showing the 

variation in Hg, Rj|, and the significance 
of the paramagnetic component compared to the 
ferromagnetic component (see Figure 11) . 

Figure 9 - Hysteresis ratios for basalts and synthetic 
titanomagnetite . 

A. Rj vs. R^ for basaltic rocks of constant 
composition but variable size modes. 

B. Rt vs. Ru for synthetic titanomagnetite 
of variable composition, but constant 
grain size. 

Figure 10 - Schematic distributions of spheres and needles 
and the corresponding magnetization curves, b 
based on the shape field for discrete compo- 
nents and experimental observations. 

A. Four possible needls-sphere configurations. 

B. Magnetization curves for each of the above 
configurations . 

Figure 11 - Schematic analysis of two components A and B 
with different shape factors. The component 
curves are normalized for simplicity. 



FIGURE I 


SIMPLIFIED MAGNETIZATION CURVES 
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FIGURE 2 


QUENCHED OCEAN BASALT 
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FIGURE 3 
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MAGNETIZATION 


FIGURE 4 
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FIGURE 5 


TWO SPHERE CONFIGURATION 

ILLUSTRATION OF THE EFFECT OF INTERACTION ON 
THE EFFECTIVE SHAPE FIELDS 
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FIGURE 7 



FUSION CRUST — FERROMAGNETIC COMPONENT IS 

MAINLY (FeNi) 3 0 4 

METEORITE BODY-MAINLY FeNi, OXIDIZED FeNi, Fe.^S 
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HYPOTHETICAL MAGNETIZATION CURVES 



FIGURE 10 





TWO COMPONENTS A, AND B OF EQUAL MAGNETIZATION 
BUT DIFFERENT SHAPE FACTORS 




REFERENCES 


Chao, E. C. T . , 1968 , Pressure and temperature histories 
of impact metamorphosed rocks-based on petrographic ob- 
servations: in Shock Metamorphism, of Natural Materials , 

edited by B. M. French and N. M. Short, Mono Book Corp., 
Baltimore. 

El Goresy, A., 1968 , The opaque minerals in impactites 
glasses: in Shock Me tamorphi sm of Natural Materials , 

edited by B. M. French and N. M. Short, Mono Book Corp., 
Baltimore. 

French, B. M., and N. M. Short, 1968 , Shock Metamorphism 
of Natural Materials , Mono Book Corp., Baltimore. 

Kneller, E. F., and F. E. luborsky, 1963 , Particle size 
dependence of coercivity and remanence of single domain 
particles: Jour. Appl. Phys., V34, p656. 

Meikeljohn, W. H., 1953, Experimental study of the co- 
ercive force of fine particles: Rev. Mod. Phys., V25, 

p302. 

Rhodes, P.,G. Rowlands, and D. R. Birchall, 1962 , Mag- 
netization energies and distributions in ferromagnetics: 
Jour. Phys. Soc. Japan, V17 (Bl), p543. 

Stacey, F. D. , J. F. Lovering, and L. G. Parry, 1961 , 
Thermomagnetic properties, natural magnetic moments, 
and magnetic anisotropy of some chondritic meteorites: 
Jour. Geophys. Res., V66, pl523. 

Taylor, G. J., and D. Heymann, 1969 , Shock reheating and 
gas retention ages of chondrites: Earth & Planet. Sci. 

Lett., V7, pl51. 

Wasilewski, P., 1972 , Magnetic hysteresis in natural 
materials: to be published in Earth & Planet. Sci. 

Lett. 

Wasilewski, P., 1972 ; Magnetic properties of chondrite 
meteorite fusion crust: NASA-Goddard Space Flight Center, 

X Document, X 644-72-163. 

Wasilewski, P., 1972 , Magnetic hysteresis in chondrite 
meteorites: NASA-Goddard Space Flight Center, X Docu- 

ment, X 644-72-167. 




nasa-gsfc 



